arXiv, 26 February 2013 



O 

(N 



Texas 78712 



Amplification of non-Markovian decay due to bound state ab- 
sorption into continuum 

Savannah Garmon 1 *, Tomio Petrosky 2 , Lena Simine 1 , and Dvira Segal 1 

1 Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, 
Toronto, Ontario M5S 3H6, Canada 

2 Center for Complex Quantum Systems, University of Texas at Austin, 1 University Station, C1609, Austin, 

(N ' 
£ 

Key words non-Markovian decay, power law decay, bound state 

It is known that quantum systems yield non-exponential (power law) decay on long time scales, associated 
with continuum threshold effects contributing to the survival probability for a prepared initial state. For 
an open quantum system consisting of a discrete state coupled to continuum, we study the case in which 
a discrete bound state of the full Hamiltonian approaches the energy continuum as the system parameters 
are varied. We find in this case that at least two regions exist yielding qualitatively different power law 
decay behaviors; we term these the long time 'near zone' and long time 'far zone.' In the near zone the 
survival probability falls off according to a t -1 power law, and in the far zone it falls off as t~ 3 . We 
show that the timescale Tq separating these two regions is inversely related to the gap between the discrete 
bound state energy and the continuum threshold. In the case that the bound state is absorbed into the 
continuum and vanishes, then the time scale Tq diverges and the survival probability follows the t _1 power 
law even on asymptotic scales. Conversely, one could study the case of an anti-bound state approaching 
the threshold before being ejected from the continuum to form a bound state. Again the t^ 1 power law 
(f) ' dominates precisely at the point of ejection. 
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1 Introduction 



For decades it has been known that deviations from exponential decay in quantum systems exist under 
quite general conditions, both on extremely short and very long time scales UJ [5J 2) ■ The short time 
scale deviations typically give rise to parabolic decay; rapidly repeated measurements of the system on this 
timescale may give rise to either the quantum Zeno or quantum anti-Zeno effects, resulting in decelerated 
| or accelerated decay, respectively J2]|5]|6]. However, these effects typically occur on quite short timescales, 

before giving way to the more familiar exponential decay, and hence experimental study was elusive for 
several decades; it is only relatively recently that these effects were experimentally verified |]7]|8]. 

Meanwhile, it was Khalfin (9), relying on the Paley- Wiener theorem ifTOl . who showed that a quantum 
system that is bounded below necessarily gives rise to non-exponential decay on asymptotic timescales; an 
alternative proof that does not rely on the Paley-Wiener theorem was provided by Hack ITU . However, the 
timescales on which the non-exponential effects manifest typically exceed several lifetimes, after which 
the survival probability is usually quite depleted, again posing an experimental challenge. It is therefore 
quite recently this effect has been verified experimentally in the luminescence decay properties of dissolved 
organic materials following laser-excitation lfl2l . The recent advances in the fabrication of nano-scale de- 
vices and cold atom trapping provide an optimistic outlook for deeper studies in this direction, particularly 
given the previous success in detecting non-exponential decay at short times for sodium atoms initially 
trapped in an accelerating optical potential Q. 

The long-time deviation, which typically takes the form of an inverse power law, is intimately connected 
with the properties of the continuum and particularly the threshold. As an example, for a particle initially 
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trapped in a shallow inner potential well it has been shown that as one varies the strength of a repulsive 
outer inverse square potential, the exponent of the power law varies in a predictable manner iTPJl . In 
another study, it has been shown for a radial potential that for the s-wave component of the wave function 
as one brings the energy of an initially prepared unstable state (a resonant state initially localized in the 
inner well) close to the threshold that the exponential decay can be completely suppressed, resulting in 
power law decay on all time scales lfl4l . A similar effect is demonstrated in Ref. 1151 as a resonant state 
approaches the threshold in a semiconductor double barrier quantum structure. 

None of the three above studies, however, consider the effect of bound states on the power law decay. 
A semi-infinite tight-binding chain with an endpoint impurity is considered in Ref. Ifl6l . in which case 
exponential decay is again suppressed, however, in this case the effect is associated with the appearance of 
a bound state. The same effect is observed in an equivalent model in Ref. [17] , as well as in a semi-infinite 
chain with a two-level endpoint impurity in Ref. ifTHl . In the latter two studies, an intensification of the 
power law decay is associated with the appearance or dissolution of a bound state (i.e. a transition from a 
localized state to a virtual state in the language of Ref. lTT8l ). We also note that the effect of a bound state 
near threshold on the non-exponential decay is considered from a more mathematical perspective in Ref. 
lfl9l . However, the analysis presented in that work is not rooted in any specific physical Hamiltonian but 
rather the authors place conditions on the spectrum itself. For example, they assume that an unperturbed 
bound state appearing directly at threshold cannot result in a bound state shifted below the threshold under 
the influence of the interaction with the continuum. However, this assumption is contradicted by both of 
the physical models considered in this paper, as well those studied in Refs. |[T6l [171 [T8ll . 

In the following study, we consider more closely the role that the bound state plays in relation to the 
non-exponential decay in open quantum systems. In particular, we show that in the case that a bound state 
approaches the energy continuum, there are two different temporal regions in which the leading term in the 
power law expansion takes on different exponents; we term these the long time "near zone" and the long 
time "far zone." Q The power law decay in the near zone is enhanced relative to that in the far zone and the 
time scale Tq that separates the two zones is inversely related to the gap between the bound state energy 
and the continuum threshold. Hence as the bound state approaches the threshold the time scale Tq grows 
larger, resulting in the amplified near zone decay extending over longer and longer time scales. Finally, 
in the case that bound state is absorbed into the continuum at the threshold (meaning that the bound state 
is destabilized as it vanishes from the diagonalized Hamiltonian spectrum) the time scale Tq diverges and 
the near zone description becomes the fully asymptotic large time behavior. 

On the other side of this transition, an anti-bound state forms in the second Riemann sheet, meanwhile 
the long time far zone behavior is restored in the asymptotic limit. Hence, one could just as well view this 
transition in reverse: an anti-bound state approaches the energy continuum, and the power law decay is 
amplified at precisely the point that the anti-bound state touches the continuum edge before being ejected 
from the continuum as a bound state. Indeed, the power law decay is easier to study in the latter case as 
neither a bound nor a resonance appear in the spectrum (at least, for the model that we consider below), 
leaving the power law as the most pronounced effect in the time evolution lfl6l [171 [181 . 

In Sec. |2] below we present our argument that the power law decay is amplified as the bound state 
approaches the continuum threshold (or the edge of the conduction band in solid state terms) and in Sec. 
I2.3l we discuss two cases for open quantum systems: in one case the bound state is prevented from reaching 
the continuum threshold due to a divergent van Hove singularity in the density of states while in the other 
the form of the interaction potential modifies the system properties and allows the bound state to reach the 
continuum edge, resulting in a divergence for the timescale Tq . In Sec. [3]we look at two model realizations 
of these effects: in Sec. 13.11 we study an infinitely long quantum dot superlattice with a side-coupled 
impurity (Model I); in this case the van Hove singularity prevents the bound state from destabilizing at the 
continuum edge. Meanwhile in Sec. I3.2l we present the case of a similar model consisting of a semi-infinite 

This naming convention is chosen in analogy with the spatial dependence in the description of the photon dressing for a two- 
level atom; in that case there is the near zone or Van der Waals zone in which the photon cloud falls off as x~ 6 and the far zone 
or Casimir-Polder zone in which it falls off as x~ 7 . Finally there is a characteristic length separating these two zones, fixed by the 
energy gap between the two levels 1201 . 
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chain with an endpoint impurity (Model II); here the bound state may reach the continuum edge, resulting 
in the divergence of the timescale Tq. We make our concluding remarks in Sec. |4] 



2 Bound state de-stabilization and non-Markovian decay 

Here we consider an open quantum system with a generic Hamiltonian of the form 

H = e d d)d + J dke k c[c k +g J dk{v k c\d + vld)c k ) (1) 

consisting of a discrete state with creation operator eft coupled to the field operator c! through the coupling 
potential gv k . The discrete state has energy e d while the continuum is given by the dispersion e k , which 
we suppose extends from a lower threshold e t ~jj to an upper threshold e^. Typically, in atomic, molecular, 
and optical systems we will only have a lower threshold while is taken to infinity, while in condensed 
matter systems both thresholds are finite as they delimit the edges of the conduction band. 

For our purposes we tend to focus in this paper on studies of the time evolution when the discrete 
energy e d is in the vicinity of the lower threshold ejjj . In this case, the upper threshold will generally have 
negligible impact on the survival probability; hence, for our purposes in this section we might as well 
simply refer to e t h = el as 'the system threshold' and assume that ejjj is largely inconsequential (or simply 
taken to infinity — > oo). 



2.1 Generic model: discrete spectrum 

We solve for the discrete spectrum of the full Hamiltonian by calculating the Green's function for the 
discrete state 

G d , d (z) = {d\^—\d) = L_^ = _L (2) 

z — H z — e d — ±j(z) ij(z) 

Here £(z) is the self-energy function that determines how the perturbed system eigenvalues are shifted 
from the unperturbed energy ed- More precisely we may write £(z) as 

E(«) = A(z)+A{z) (3) 

in which A(z) is analytic for all z, meanwhile A(z) represents the non-analytic component of £(z); specif- 
ically A(z) is undefined along the branch cut associated with the continuum e k . Instead, one must define 
this function on the two different Riemann sheets associated with the branch cut as A I,u (z) = A(z ± i0 + ). 
Then the bound state solutions zsi of the system, for example, are written as solutions of 



"-B. 



e d - E 1 ^) = (4) 



in the first sheet. 

Further, we have defined in Eq. (O the dispersion function ?y(z) = z — td — E(z), which yields the 
discrete spectrum as the solutions of the equation r\{z) = 0. For many simple systems we may explicitly 
obtain rj(z), which then allows us to write the dispersion relation in polynomial form with the eigenvalues 
given as the roots of the nth order polynomial equation 

P n (zi) = SfeWizdrpfe) - 0, (5) 

with the n resulting eigenvalues Zi(ed,g) functions of the Hamiltonian parameters e d and g; for this paper 
we will usually assume we may vary e d while g takes a fixed value and may be repressed in the notation 
from this point forward. Note that S(z) in Eq. (0 is simply a z-dependent factor that might need to be 
multiplied through to obtain a polynomial. Let us further assume that for values €d < (a there exists 
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Fig. 1 Complex z plane with the branch cut extending from the threshold e t h to oo along the real axis. A single 
bound state zb appears below the threshold. The Green line represents the (a) initial contour (Bromwich path), and (b) 
deformed contour, with unbroken lines in the first sheet and dashed lines in the second sheet. 



a single bound state in the spectrum with eigenvalue z B (e d ), which satisfies z B (e d ) < e th (since we are 
dealing with an open system, the other eigenvalues might be any combination of anti-bound states, resonant 
states or anti-resonant states). However, if we increase the value of e d such that e d — > £a then consequently 
the bound state energy approaches the threshold zg — > eth. At the precise value e d = £a the bound state 
reaches the continuum threshold zb = e tI1 and is absorbed H211 . Hence the bound state is no longer present 
in the system for values e d > e A (instead it becomes an anti-bound state for e d > e A ). Our resulting 
diagonalized Hamiltonian hence takes the form l22l 

{z-e,d\d B + / dke k c\c k e d < £a 

(6) 
/ dke k c\c k e d > e A 

Note that for our purposes in this paper, when we talk about 'destabilization of the bound state,' we are 
referring to the transition from the first line to the second line of this equation^ Finally, for ease in the 
following calculation, we re-write the bound state eigenvalue in terms of Aq, the energy gap between the 
bound state eigenvalue and the continuum threshold, as 

ZB(e d ) = ea, - Aq. (7) 



2.2 Generic model: survival probability 
Now we calculate the survival probability P(t) = \A(t)\ 2 in which 

A(t) = (d\e- iHt \d) = / dze- izt {d\-^—\d) = -L / dze-«*-l- (8) 
Zirt J T z — ti im J T ly [zj 

is the survival amplitude for the system in the prepared initial state \fy) = \d). Meanwhile, r] 1 (z) is the 
dispersion function defined on the first Riemann sheet and which may be analytically continued into the 
second sheet on the lower half of the complex energy plane through the branch cut associated with the 
continuum e k . Finally, the contour T is the Bromwich path 12311241 that extends just above the real axis in 
the first sheet as [— oo + iS, oo + iS] as show in Fig. [TJa). 

As a practical means for evaluating the survival amplitude we now deform the contour by the standard 
approach, which is to drag the entire contour out to infinity in the lower half plane. In the case e d < £a 



Also note that other transitions might occur for values of ed much larger or much smaller than eA (appearance of resonance 
state, etc.) however we primarily focus on the vicinity of the transition ~ £a in this paper. 
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then we will pick up a pole associated with the discrete energy zb of the bound state located on the real axis 
below the threshold e t h. The remaining contour consists of a portion at infinity that gives no contribution 
and a segment T\ that runs from the threshold e t h out to infinity at e t h — ioo. Meanwhile the segment of the 
original contour that is initially located just above the continuum must be analytically continued into the 
second sheet through the branch cut. As we pull the contour down, we generally pick up a pole for each 
resonance state that might be present, as well as a second contour Tu that extends from e t h out to infinity 
at eth — ?'oo in the second sheet. Hence ultimately the contributions to the survival amplitude are 

A(t) = A B (t) + A ± (t) + A R (t), (9) 

in which A B (t) = Rse^ 1ZBt is the contribution from the single bound state with residue Rr, A^(t) is the 
sum of contributions of any resonance poles (for simplicity in our example, let us assume there are none) 
and then the contribution from the two contours T\ and Tn take the form of the so-called 'background 
integral' 

4h(*) = t — / dze~ lzt ^—^ + — / dze- lzt -r-- 

= i r~ to ° dze -^A^i S ( Z ). do) 

WJe* P n{z) 

Now, let us explicitly write the polynomial P n (z) in terms of the single bound state as 

P n {z) = {Z - Z B ) Rn-l(z) (11) 

with R n -i(z) an n — 1 order polynomial. Performing a change of the integration variable according to 
z = £th — iy we may write the integral Eq. ([Tot as 

f°° rin p—yt- A (^-w) 



A th (t) = -Zie"*** / dye-yt— ^ — . (12) 

From this point, we make a final coordinate transformation s = yt and re-write the polynomial in the 
denominator as 

2 

— s s s s 

i? n _i(-) = i?„_i(e th - i-) = f +fi- + f 2 '-pr + ■■■ (13) 
t t t t z 

Finally we obtain the background integral in the standard form as 



2e -i*t f°° j _ _ s A( et h-if) 



it Jo A Q r + Zif +l 2 jr + ... 

in terms of a new polynomial L n (s/t) = Iq + sli/t + . . . with lo = AqTo- 

Having obtained Eq. ( fT4b . we are now in a position to make some comments about the non-Markovian 
decay. First, note that for small t, the rightmost terms in the denominator of this integral will make the 
greatest contribution to the integration. However, we are interested in the long time effects; let us first 
focus on the temporal region 

£<t<-^ (15) 

When this condition is satisfied, the second term in the denominator of Eq. ( TBI will dominate the time 
evolution. We may expand the denominator accordingly in order to write the approximation 

4h(t) * A^(t) = — / F(t'; t)dr' (16) 
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in which we have defined the function 

F(r;t)= / ds e~ ST A{e th - i-), (17) 
Jo 1 

which will appear again momentarily. The label "NZ" in the above equation stands for long time near 
zone. 

Now let us consider the fully asymptotic case 

-rV « t. (18) 

Here the first term in the denominator of Eq. (TBI will dominate the integration, resulting in a new approx- 
imation in the form 

4h(i) « A™(t) = -^—F(r = l-t). (19) 

for the long time far zone. 

Now let us compare the approximations in the two different zones in Eqs. ( fTol l and (fl~9b . Notice that if 
we assume 



F(T';t)dT' ~ lim F(t = l;t). 



(20) 



then the two zones obey differing powerlaws A™(t) ~ j^4|J z (i), which in terms of the survival probability 
itself yields 

(*»!)• (2D 

Hence we find a relatively amplified power law decay in the near zone in comparison to the far zone. 
However, if we consider the case that Aq = such that the bound state is absorbed into the continuum, 
then the time scale Tq = Ii/AqTo in Eq. ( fT~8b diverges, and the relatively amplified near zone decay 
becomes the full asymptotic behavior of the system (that this is true is fairly easy to see from Eq. (TBI as 
well). 



2.3 Comment on bound state de-stablization in OQS 

In this section we comment on the behavior of bound states in open quantum systems and the de-stabilization 
process itself. In the previous section we realized a process by which amplification of the long time power 
law decay may occur in such systems as the bound state is absorbed into the continuum. Here let us remark 
on the circumstances under which such a process occurs. 

The key point is that this behavior is largely determined by the form of the non-analytic portion of the 
self-energy A(z) from Eq. ||3}, which contains the branch cut associated with the continuum. Hence in the 
single band case when the bound state is absorbed at the edge of the continuum we have A(zb = £th) = 0. 
Now the form of A(z) in turn is largely fixed by two pieces from the Hamiltonian: (A) the interaction 
potential i>k, and (B) the continuum density of states p(ej,) (proportional in the 1-D case to the derivative 
of the inverted dispersion relation dk/dek)- As shown in Ref. l29l , in the case 

lim K( e )| 2 p(e) -> oo, (22) 

then the condition A(e t h) = can never be satisfied in Eq. ©, which prevents the bound state from ever 
'touching' the continuum and hence de-stabilization never occurs. Instead, as ^> e t h then the bound state 
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will become 'stuck' on the outer edge of the continuum, and will only be able to approach the threshold 
asymptotically (we refer to these type of bound states as persistent bound states as in Ref. |29l ). 
On the other hand, if we consider a system in which 

[Ke)| 2 p( e )U, h =0 (23) 

is satisfied, then the condition A(e t h) = may be satisfied for some choice of system parameters at which 
the bound state will destabilize through absorption into the continuum, ultimately producing an anti-bound 
state in the second Riemann sheet (or conversely, an anti-bound state may approach the continuum before 
being ejected as a bound state in the first sheet) ||2T1 . This is equivalent to the localized to virtual state 
transition studied in Ref. 11251 . 

Below we will study two models that provide examples of these two cases. In Sec. I3.1l we will consider 
an infinite chain with a side-coupled impurity, which satisfies the relation in Eq. d22l . As a result we find 
that the persistent bound state may move closer and closer to the band edge, but the condition Aq = 
cannot be satisfied for finite values of the system parameters. Hence the time scale in Eq. ( TT~8T > never 
approaches infinity and the long time far zone behavior can never be entirely eliminated from the system. 

However, in Sec. I3.2l we study a semi-infinite chain with an endpoint impurity, which satisfies instead 
the relation Eq. d23l . In this case, it is possible to vary the system parameters such that Aq = and then 
the amplified power law decay associated with the long time near zone will dominate the long time system 
dynamics even in the asymptotic limit. As a further point, both bound states and resonance states may be 
eliminated from the system entirely in this model for a wide range of parameter values (that includes the 
special point Aq = 0), resulting in non-exponential evolution on all time scales. 



3 Model realizations of amplified non-Markovian decay 

We now present realizations of the non-Markovian amplification effect relying on two variations of tight- 
binding chains with coupled impurity sites. Our two models may be compactly written as 

H = e d cPd - t ^(4c„+i + 4 +1 c„) + e f _ l ^2 c n c n ~ 9 (4 d + ^ c n ) (24) 

n n 

in which c„ and d are the annihilation operators along the chain and at the impurity site, respectively, while 
t sets the tunneling strength along the chain and g gives the coupling strength to the impurity. Meanwhile 
e M gives the chemical potential at each site along the chain and is the impurity site energy. We fix the 
energy scale according to units t = 1 and set the origin as = 0. Note that the creation operator is 
the creation operator for the site coupled to the impurity (uq will be specified below). 

For Model I, we have the infinite geometry chain with N + 1 sites, indexed as n = — -y , . . . , ^ , with 
periodic boundary conditions and an impurity side-coupled at the middle of the chain no = as shown in 
Fig. 12 a). For Model II, we have the semi-infinite geometry chain with N sites, indexed as n = 1, . . . , ^, 
with open boundary conditions and an end-point impurity attached at the no = 1 site as shown in Fig. [2jb). 

We visualize that these models may be physically realized as either a semiconductor superlattice as in 
Ref. 1261 or a series of coupled single mode waveguides as in Ref. lfl6l . 



3.1 Model I: persistent bound state supported at threshold by DOS divergence 

Let us now consider Model I, the infinite-chain realization of the Hamiltonian in Eq. ([24). After performing 
a Fourier transform on the chain sites c n and taking the limit iV^oowe may write 

H = e d dU+- dke k c f k c k +g v k dk (c\d + $ c k ) (25) 
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Fig. 2 (a) Model I: infinite geometry with side-coupled impurity, and (b) Model II: semi-infinite geometry with 
endpoint impurity. 



in which the dispersion for the diagonalized modes is given by = — cos k and the potential Vk = 
— 1/ V^r is simply a constant. We note immediately the form of the normalized density of states for our 
dispersion takes the form 

, N 1 1 dk 

p(tk) = , 5 oc j—- (26) 

with divergences at either band edge = ±1. Hence at either band edge the condition 

lim , \ v k(e)\ 2 p(t) -> oo (27) 

is satisfied, which is the equivalent of the condition previously appearing in Eq. (1221 . As a result, the 
bound state absorption will never be realized in this system, and the amplification of the power law decay 
will only occur up to a certain timescale. 

3.1.1 Model I: spectrum 

The discrete spectrum of Model I has been previously studied in Refs. [25, 26, 27, 28, 29). As presented 
in Ref. l26l . the spectrum may be obtained as the poles of the Green's function for the impurity site, which 
after re-summation appears as 

G d , d {z) = {d\-^—,\d) = 1 ( 28) 

z-H z-e d -T>(z) r)(z) 

with the self-energy function S(z) given by 

Z(z)=g 2 f-^dk=^=. (29) 
J z-e k ^Jz 2 - 1 



Note there is only one term in the self-energy, such that A(e) = g 2 / 'y 1 z 2 — 1 and A(z) = 0. Further we 
see immediately that A(±l) = can never be satisfied, as a direct result of Eq. (127V 

The discrete spectrum is now fixed by the poles of the Green's function as the solutions to rj(z) = or 

g 2 

z-e d = ^==. (30) 

By squaring Eq. (130b we immediately obtain an equivalent quadratic polynomial equation for the spectrum. 
We plot the resulting real part of the spectrum as a function of e d in Fig. EJa), in which solid lines represent 
the two bound state solutions appearing in the first Riemann sheet and dashed lines give the two other 
solutions appearing in the second sheet. 

For our purposes, let us focus on the region of the spectrum in the lower half of Fig. |3ja), in which the 
impurity energy lies fairly close to the lower band edge e d ~ — 1. We show this region in greater detail 
in Fig. |3b). For the region in which e d < — e 7 (the approximate value of e 7 is given momentarily), the 
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Fig. 3 (a) Discrete spectrum for Model I as a function of td with g = 0.4, and (b) close-up in the vicinity of the lower 
band edge. The two bound state solutions z± are both shown as full, unbroken (blue) curves while anti-bound states 
and resonance states in the second Riemann sheet are shown with dotted curves. The lower and upper band edges at 
y = ±1 as well as the unperturbed impurity energy y = td are shown with unbroken (black) lines. The values for 
±e 7 where the resonance state pair forms are also indicated (see Eq. <3U for approximate value). 



spectrum consists of two bound states, one appearing in the lower half of the spectrum with discrete energy 
z_ ~ e d , and a second on the positive half of the spectrum lying close to the upper band edge z + ~ 1, 
as well as two anti-bound states in the second Riemann sheet. When we calculate the survival probability 
here, not surprisingly we find that the lower bound state z_ dominates the time evolution while the upper 
bound state has a largely negligible effect. 

If we choose to gradually increase the value of the parameter e^, we find that eventually two new features 
occur in the spectrum. First, at a specific value of the impurity energy €d = — e 7 , indicated in Fig [3] the 
two anti-bound states collide on the second sheet to form a resonance, anti-resonance paii0 with complex 
conjugate eigenvalues. This complex conjugate pair persists for all values |e<j| < e 7 with the approximate 
value 

e 7 = l + §ff 4/8 + 0(ff 8/8 ) (3D 

obtained in Ref. 11261 (a plot of the decay width may also be found in that paper). Second, as we continue 
to increase e c [ such that it exceeds values on the order of > —1 — g 4 ^ 3 , the lower bound state takes 
on a value Z- < —1 as it gets 'stuck' on the lower band edge, forming a persistent bound state. Q That 
this happens can be seen intuitively from Eq. J30t ; if z 2 — 1 ~ g 4 is satisfied in the root on the RHS, we 
see that both sides of the equation may become order unity. More precisely, by following the perturbative 
calculation detailed in Ref. ||29ll we obtain 

= -l-A Q +0(g 8 ), (32) 

for all values ed S> — 1. Notice here that z_ never actually reaches to the band edge at —1 except in 
the limit e c i — > oo (including higher order terms would be consistent with this statement). Hence, for the 
present model, the bound state is always maintained, hovering just outside the edge of the continuum, as 
we predicted from the condition reported in Eq. (|27V We have also defined Aq here as the gap at lowest 
order between the bound state and the lower band edge. As we evaluate the survival probability below, 
we find that the background integral naturally appears written in terms of Aq, rendering unnecessary an 



^ The resonance, anti-resonance pair forms at a real- valued exceptional point I30ll31ll32l , a process which is described in greater 
detail in Ref. 1331 ; for a clever experimental demonstration of this process applying simple LCR circuits, see Ref. 1341 . 
The transition to the persistent bound state appearing at the band edge is described in detail in Ref. 1291 . 
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explicit factorization of the polynomial dispersion in terms of the bound state eigenvalue (as performed 
previously in Sec. \2.2\ . 



3.1.2 Model I: Survival probability 

We now calculate the survival probability along the lines of our previous calculation in Sec. 12.21 Following 
the steps leading from Eq. ^ to Eq. (TTOt we obtain an expression for the background integral contribution 
to the survival amplitude for Model I as 

2 ._l-ioo / 2 _ 1 

AhW = -. dz e— * - 1* j. (33) 



7riy_! (z-e d y(z 2 -l)-g* 

Next we make the transformation z = — 1 — iy in which y runs from to oo, followed by the second 
transformation s = yt in order to obtain 



y/ s 2 — 2its 

Ah(t) = "\ I dse~ s =r^r 2 ^Fa- (34) 

*g 2 t 2 Jo i - i S lz + s 21r£i + ls 3l<gk - s ilan 

Here we have defined the key time-scale 

_ 2{l + e d f _ 1 

q = ^^~Aq~ (35) 

as well as T 2 = (5 + e rf )/2(l + e d ), T 3 = s/2 + T d /{\ + e d ), and T 4 = (2(1 + e d ))~ 2/3 . Now in the case 
of T 2 < i < Tq, we have 

_, Vs 2 - 2iis 

tt^Tq Jo 



git 



F (t', i) cIt' 



irg 2 tT Q J T=1 
2 it 

*(-l/2,0;-2ii) (36) 



20F(l + £d) 2 t 

in which ^(a, c; z) is the confluent hypergeometric function and 



— , . f°° , ST ^ — - — ite- ltT KA-itr) 

F(t; t)= ds e' s Vs 2 - 2iis = ^ '- 

Jo T 



— itr i 

■ 

(37) 



with K\ (x) the modified Bessel function. As we generally expect the condition T2 <C t *C Tq will apply 
for t> lwe apply the asymptotic expansion for the confluent hypergeometric function ('5 (a, c; z) — > z~ a 
for z — > oo) to obtain 

„37ri/4 2 it 

< z (i) « ° , g ° . (38) 



^F(l + e d ) 2 ^/2 

Hence we find the survival probability goes like P NZ (i) ~ t" 1 in the long time near zone. Meanwhile 
from Eq. ([34b in the long time far zone t 3> Tq we utilize the asymptotic behavior K\{x) — > yj-n /2xe~ x 
to obtain 



4h Z W « — ^5^^*) 
ng z t z 



—ifi(-it), (39) 
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Fig. 4 Numerical results for the time evolution P(t) for Model I for the three cases. 



which in the asymptotic limit gives 

Af(t) » 1=4— . (40) 



2ngH 3 / 2 

Using Eqs. ([38) and (00) to calculate the survival probabilities P^ z = — t" 1 and P^ z = 

IAh Z (0l ~ * _3 > we immediately verify that the predicted relation Eq. (EH holds true for the present model. 
Unfortunately however, the fact that both bound states and a resonant state are present in the system for the 
parameter choices we are interested in (i.e. Aq small) results in a decay profile in which it is difficult to 
clearly distinguish the time evolution behavior in the near zone; this is demonstrated by numerical results in 
Fig. [4] For the case g = 0.1, = —0.5 (large red diamonds) we see that the exponential decay dominates 
the visible time evolution up to about t < 20, then the bound state is the dominant feature around t > 50. 
In the approximate region 20 < t < 50 the evolution evolves according to interference between these 
two terms, resulting in an oscillation pattern. Then for g = 0.2, = —0.5 (blue dots) the decay width 
(proportional to g 2 ) is smaller resulting in an elongated timescale t < 60 for the exponential decay. Finally 
for g = 0.1, €d = —0.5 the decay width is even smaller, resulting in an exponential decay that outlives the 
plot parameters. 

We see that the power law decay effects are not apparent in any of these plots. Fortunately, we may 
overcome this difficulty for certain parameter choices in the case of Model II. 

3.2 Model II: bound state absorbed into continuum 

We now turn to the Model II realization of the Hamiltonian Eq. (124) . consisting of a semi-infinite-chain 
with an endpoint impurity. We perform a half-range Fourier transform on the chain sites c n and take the 
limit N — > oo in order to write 

H = e d dU+- dke k c\c k +g v k dk (c\d + ef 1 c k ) (41) 
^ Jo Jo 

with the same tight-binding dispersion e k = — cos k as for Model I, hence the density of states takes 
precisely the same form as Eq. (|26) . including the square root form of the divergences at either band edge 

e k = ±1 

Meanwhile the potential takes the form v k = — sink/y^n, which introduces some quite different prop- 
erties into the system in comparison to Model I. The primary point is that at either band edge the condition 
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Fig. 5 (a) Discrete spectrum for Model II as a function of ed with g — 1/2 (linear case). The lower and upper band 
edges as well the unperturbed energy y — ed are all shown with solid black lines. The bound state in the first sheet 
is shown by a full, unbroken (blue) curve while the anti-bound state in the second sheet is represented with a (blue) 
dotted curve. Vertical dashed lines represent the point = ±1/2 at which the bound state brushes the continuum 
to form an anti-bound state (or vice versa), (b) Numerical results for the time evolution P(t) for Model II for two 
cases: ed = —0.4 (large red diamonds) along with the red line representing the far zone evolution P FZ (t) ~ i -3 , and 
ed = —0.5 (small blue dots) along with the blue line representing the near zone evolution P NZ (t) ~ t -1 . 



[k(e)| 2 p(e)] e=±1 = (42) 

is satisfied, which is the equivalent of the condition Eq. d23l that allows for a bound state to become 
destabilized as it touches either edge of the conduction band at specific parameter values. As we study the 
spectrum in greater detail in the next section we will also present the precise value = €a where this 
destabilization occurs (assuming a fixed value of g). 



3.2.1 Model II: spectrum 

The discrete spectrum and the power law decay for Model II has been studied previously in Ref. Ifl6l for 
the case of e c i — (all sites with equal energy) while only the coupling to the endpoint may be varied. For 
our purposes of illustrating a simple generic feature we will focus instead on the particular case of a fixed 
endpoint coupling (g = 1/2 in Fig. 12(b)) while the impurity energy ed may be varied. The power law 
decay has been studied in models similar to these in Refs. ifTTl [T8l l24l . Just as for Model I, we obtain the 
discrete spectrum for the present model from the Green's function at the endpoint impurity site as in Eq. 
(1281 1. And as in Eq. (29) , it is a straight-forward exercise to obtain the self-energy function for arbitrary g 
as 

E(z) = 2g 2 (z - Vz 2 - l) (43) 

for the semi-infinite chain model, in which it is now possible for A(±l) = to be satisfied for A(z) = 
2g 2 Vz 2 ~T. 

The discrete spectrum is obtained from 77(2) = z — — S(z) = as 

z-e d =^(z- Vz 2 - l) (44) 

for the specific case g = 1/2. In the more general case of arbitrary g, we would obtain an equivalent 
quadratic polynomial after squaring out the root on the RHS of this equation. However, for our present 
case we instead find a simple linear equation. Before presenting the solution, we first note there is a special 
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point in parameter space associated with the values z = ±1 at which the root in Eq. (l44l vanishes. What 
happens at this special point is that a bound state appearing in the first Riemann sheet touches the edge 
of the continuum, before giving rise to an anti -bound state in the second sheet (or vice- versa) ifTH] 1211 . 
Plugging the values z = ±1 into Eq. (f44b yields e d = ±1/2 as the precise points in parameter space where 
these transitions occur. It can be further demonstrated that for \e d \ > 1/2 there is exactly one bound state 
present in the system, and for \e d \ < 1/2 there is exactly one anti -bound state. In either case, the linear 
dispersion yields the single solution as 

z L (e d ) =e d +—, (45) 

which holds for the bound and anti-bound state alike. We plot the spectrum as a function of e d in the Fig. 
|5ja), in which unbroken curves show the bound state portion of the spectrum, while the dotted curves show 
the anti -bound state. It's interesting to note that in the anti -bound case le^l < 1/2, there is neither a bound 
state nor a resonance state in the system, hence the time evolution is non-exponential on all time scales 

MM- 



3.2.2 Model II: Survival probability 

Again we write the survival probability for Model II by following the standard method from Sec. 12 . 2 1 to 
obtain 



— 1 — 200 



A th (t) = - 1 / dz e 

2nie d 



l 





- 1 


z — 




Vs 2 - 


2its 



dse s =r— , (46) 



2me d A Q P 7 1 + is^ 



following the usual transformations z = — 1 — iy then s = yt. Here Aq = z^ + 1 is the gap between the 
bound state and the lower band edge at z = — 1 and Tq = (Aq) -1 is the usual time scale dividing the 
near zone and far zone. In the long time near zone 1 t <IC Tq we have 

■it i-oo 

* -l^r- / F(r',t)dr' 
2me d t J T=l 

„it 

-*(-l/2,0;-2^) 



2 v / 7re ( zi 

it 

(47) 



37ri/4 g it 



2 7 re d t 1 / 2 

with the usual t -1 / 2 behavior. Meanwhile as we turn to the long time far zone in the fully asymptotic case 
t > Tq, we find from Eq. gUl that 



7ri/4git 



2V2^e d A Q t 3 / 2 ' 



(48) 



Clearly we see that in the case that the bound state touches the continuum at e d = 1/2 (Aq = 0), then 
Tq —> oo and the long time near zone gives the dominant behavior even in the asymptotic limit. 

We illustrate this by numerical results for P(t) in Fig. EJb); in the case that e d = —0.4 (red diamonds), 
there is no bound state present in the system and we clearly observe the transition from the near zone with 
slope coefficient ~ — 1 (blue line) and the far zone with coefficient ~ —3 (red line). Meanwhile for the 
value e d = —0.5 (small blue dots) we have realized the case Aq = 0, hence the far zone time scale has 
been pushed out and the near zone behavior dominates the evolution even on larger time scales. 
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4 Conclusion 

In this paper we have demonstrated in simple terms that there is some connection between the proximity of 
a bound state to threshold and the precise details of the non-exponential decay of the survival probability at 
long times in open quantum systems. This connection is encapsulated in Eq. (fT~4T > in which it is clear that 
closing the gap Aq between the threshold and bound state should have some effect on the long time power 
law decay (most likely, it will be enhanced). Further the existence of the gap gives rise to a time scale that 
divides the long time system dynamics into two zones: the long time near zone and the long time far zone. 
The question of how the power law decay will compare in these two zones largely depends on whether the 
condition in Eq. ( T20t holds. 

In Sec. [3] we explored this question for two specific nearest neighbor models with attached impurity 
levels. In these two cases Eq. ( f20b indeed held, and hence the power law decay was enhanced in the 
long time near zone in comparison to the far zone according to P FZ ~ -pP NZ as predicted in Eq. (l2TT i. 
Further, for Model I the divergent van Hove singularity in the density of states at the edge of the continuum 
prevented the bound state from being able to reach the band edge and closing the gap Aq (that is to say, 
there was a persistent bound state as defined in Ref. 11291 ). Hence one could only approach the regime 
in which the near zone dominates the power law decay for very large times only in the limit eg — > 00, in 
which case the decay contribution in Eq. (l38l would vanish anyways. However, for Model II the precise 
form of the interaction potential essentially cancels out the divergence from the density of states, which 
results in a spectrum in which it is possible to close the gap Aq = for precise values of the system 
parameters. Doing so gives rise to a situation where the enhanced power law decay associated with the 
long time near zone pushes out the far zone behavior entirely and dominates the time evolution in the fully 
asymptotic limit. 

Interestingly, it has also been demonstrated in open quantum systems that the same crossover from a 
bound state to an anti-bound state at the continuum threshold gives rise to a maximum in the local density 
of states response function at the impurity site 11251 . Indeed, this result is closely related to the amplification 
of the power law decay at threshold [ 18 1. 

However, the precise details of the transitions in the asymptotic decay behaviors for a generic system 
remains an open question. On the one hand, the t~ 3 behavior in the far zone appears to be a rather generic 
feature J2). On the other hand, one would hope to obtain a more general understanding of the near zone 
behavior. One might gain insight into this question by considering what more general statements might 
be made about the condition expressed in Eq. d20l i and under what conditions one might expect this 
relation to hold. One interesting proposal would be to study the non-exponential decay dynamics in a two- 
dimensional system, in which the inverse square root singularities ubiquitous in one-dimensional systems 
would be replaced with logarithmic divergences. This will be a future direction of study. 

Finally for completeness we note that even the typical exponential decay itself experiences anomalous 
behavior in the vicinity of the threshold. As has been pointed out in quite different physical contexts, 
Fermi's golden rule, which predicts that the decay rate to lowest order approximation is proportional to the 
density of states, breaks down in the vicinity of a singularity in the density of states 1261 1291 551 . 
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